Reconstructing the nucleon-nucleon potential by a new coupled-channel inversion method 
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A second-order supersymmetric transformation is presented, for the two-channel Schrodinger equation with 
equal thresholds. It adds a Breit-Wigner term to the mixing parameter, without modifying the eigenphase shifts, 
and modifies the potential matrix analytically. The iteration of a few such transformations allows a precise fit 
of realistic mixing parameters in terms of a Pade expansion of both the scattering matrix and the effective-range 
function. The method is applied to build an exactly-solvable potential for the neutron-proton Si-'^Di case. 
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In quantum scattering theory, the fundamental inverse prob- 
lem consists in deducing the interaction potential between two 
colliding particles from their experimental elastic-scattering 
cross sections These cross sections have first to be 

parametrized in terms of energy-dependent partial-wave phase 
shifts or scattering matrices. For a central interaction V{r), 
the partial waves decouple and a sequence of single-channel 
inverse problems have to be solved. For more complicated in- 
teractions, like the tensor interaction in nuclear physics, par- 
tial waves may be coupled and matrix potentials have to be 
constructed from coupled-channel scattering matrices. 

A formal solution to these inverse problems can be written 
in terms of integral equations [1]. In practical applications, 
experimental data turn out to be precisely parametrizable in 
terms of separable kernels for these equations. These coiTe- 
spond to scattering matrices that are rational functions of the 
wave number The integral equations can then be solved an- 
alytically and the coiTesponding potentials are also expressed 
in a separable form. This procedure applies to both the single- 
and coupled-channel cases fllltl. In the single-channel case, 
the same potentials can be more efficiently constructed with 
the help of supersymmetric quantum mechanics (SUSYQM) 
10, Eh 5 which directly relates the potentials to their scattering- 
matrix poles. Moreover, a small number of poles is in gen- 
eral sufficient to fit experimental data on the whole elastic- 
scattering energy range; this is in particular the case for the 
neutron-proton singlet (spin 0) channels , that decouple 
because of the vanishing tensor interaction. 

The present Letter aims at extending this very efficient sin- 
gle-channel method to the two-channel case without thresh- 
old difference, and at applying it to the neutron-proton triplet 
(spin 1) coupled channels. This system was studied in the 
framework of the integral-equation method by Newton and 
Fulton at low energy, and by Kohlhoff and von Geramb 
Jstl, on the whole elastic-scattering range. The present paper 
subsumes these works, by separating the effect of the coupling 
between channels in the inversion procedure, by parametriz- 
ing the data with a minimal number of poles, and by deriving 
simple expressions for the coiTesponding matrix potential. 

We consider two channels with equal thresholds and angu- 



lar momenta / and I + 2. The scattering matrix is parametrized 
by the eigenphase shifts 5i, 62 and the mixing parameter e, 

S(k) = R[e{k)] diag (^e^^s,(k)^^2^5,(k)\^ R^[e[k)] . (1) 

Here, k is the (complex-valued in general) wave number, cor- 
responding to the center-of-mass energy E = k^ in reduced 
units, and R is an orthogonal matrix 



i?(e) = 



cos e sm e 
— sin e cos e 



(2) 



We aim at building a symmetric interaction matrix V with en- 
tries Vii, V12, V22, by inversion of the scattering data 61,62, e. 
In Ref. 1 8], we show how to split this two-channel inverse 
problem into two independent parts: (i) fitting the eigen- 
phase shifts independently for each channel, with the stan- 
dard single-channel method 151-^], and (ii) fitting the mixing 
parameter, with a new type of eigenphase preserving (EPP) 
transformation. In Ref. (81], such a transformation is intro- 
duced but it is restricted to a mixing parameter which is an 
odd function of the energy, as in Ref. |2]. Here, we overcome 
this restriction, which allows us to fit more realistic mixing 
parameters. 

Let Hq = —l2(P/dr^ + Vq be an initial 2x2 matrix 
Hamiltonian, where I2 is the identity matrix and the po- 
tential Vo{r) is symmetric. With two successive SUSYQM 
transformations, this Hamiltonian can be transformed into an- 
other Hamiltonian H2, which has the same eigenphase shifts 
as Hq but a different mixing parameter Equivalently, this 
can be achieved with a two-fold EPP transformation defined 
by the intertwining relation LHq = H2L, where i is a 
second-order differential matrix operator, called transforma- 
tion operator, which maps solutions of both Hamiltonians as 
'i'2{k,r) — L'i'o{k,r). Operator L and its adjoint obey 
the factorization properties L^L = {Hq — Ei){Hq — El), 
LL^ = (H2 ~ Ei){H2 - El), where Ei ^ k^ ^ Ejn + iEn 
and El = Em — iEn are mutually conjugate factorization 
constants. General properties and different equivalent expres- 
sions for L can be found in Ref. ijstl . 
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The transformed potential reads 

V2{r)^Vo{r)~2W!,{r). 



(3) 



The matrix W2 is expressed in terms of a complex matrix 
U, made up of two vector solutions of the coupled-channel 
Schrodinger equation, -U" + VoU = EiU, as [li] 



W2{r) = -En {lMU'{r)U-\r)]] ' , 



(4) 



Potential V2 is real, symmetric, and regular when the first vec- 
tor of U is regular for r ^ Q and its asymptotic behavior reads 



r — )■ 00 — )■ 



(5) 



As announced, for this EPP transformation, the eigenphase 
shifts of the transformed scattering matrix 5*2 coincide with 
the initial ones: 6c:2{k) = Sc-.oik) for channel c = 1,2. On 
the other hand, the mixing parameter is modified as 



e2{k) = eo{k) + arctan 



Eji 



Em - k^ ■ 



(6) 



which corresponds to additional 5-matrix poles in fc^ = 
Em ± iEji. In Ref. |8], this result was proved in the case 
of purely imaginary poles, Em — 0. We made this choice 
because a mixing parameter should vanish at zero energy (up 
to an unimportant integer multiple of tt/2) for the potential 
to be short ranged (except for the r^^ centrifugal diagonal 
term). For Em 7^ 0, all calculations of Ref. |8] are actually 
valid but lead to a non-vanishing mixing parameter, £2 (0) = 
arctan(£'/i/i?fli) = ai, and to r^^ off-diagonal terms in 
the potential. Here, we solve these problems by applying an 
energy-independent rotation after the EPP transformation: the 
S matrix 5*2 transforms into S2 = {ai)S2R{ai), which 
shifts the mixing parameter by a constant value, e2{k) = 
e2{k) — Q!i, whence £2(0) = 0. Meanwhile, H2 transforms 
into H2 = R'^ {ai)H2R{ai), which has no r^^ off-diagonal 
term. 

Iterating P times such EPP transformations, with energies 
Ej = Ejij ± iEjj , j = 1, . . . ,P, one builds a chain of Hamil- 
tonians. Ho — > H2 . . . ^ H2P- The final potential V2P 
corresponds to the eigenphase equivalent S matrix 

S2p{k) = i?(e2p)diag ^e^^s,.,o ^^2^5,.,o^ R^{e2p), (7) 
with e2p{k) = e2p{k) — e2p(0) and 



f-2p{k) — eo(fc) + arctan 



E13 



Epj 



(8) 



As above, to get a vanishing zero-energy mixing parameter, 
we perform a compensation rotation of angle — e2p(0) = 
— ^^-^-^ Qfj, with ckj = arctan{Ejj/Eiij), j = 1,...,P, 
on both the scattering matrix and the Hamiltonian H2P- 

The sum in parametrization ([8]) is flexible enough to fit real- 
istic mixing parameters, even with small values of P. This can 



be explained with an effective-range expansion: the elements 
of matrix are analytical at the origin and can be Taylor ex- 
panded as a function of the energy 191]. Consequently, tan e is 
well approximated by a Pade expansion, which is exactly what 
our EPP transformations lead to: for eq = 0, tan e2P, with 2P 
arbitrary parameters, is the most general Pade approximant of 
order [P/ P] in energy, vanishing at the origin. 

This is very analog to the standard single-channel 
SUSYQM inversion |5-7], where Pade approximants also oc- 
cur There, the one-channel S matrix is expanded in terms of 
its poles iK„i in the complex wave-number plane as 

m 

Ml . , 

Si{k) = e^'^'(^-) = n {Ml >2l + 1), (9) 

m— 1 

which corresponds to the single-channel phase shift 5i{k) = 
— X]m=i arctan(fc/K„i). The poles are generally restricted 
to the imaginary axis, except for resonances, and have to be 
chosen so that the effective-range function 



Ki{k^) = cot 6i{k) = ifc2'+ 



1 



Si{k)~l 



(10) 



does not vanish at zero energy, which corresponds to a finite 
scattering length. At higher orders, Ki is usually expanded as 
a Taylor series, but this expansion often breaks down at high 
energy. The S matrix (|9]i, in contrast, leads to a Pade approx- 
imant for Ki, of order ^ — Z — l] in energy for even 
Ml and of order [^^/^^ - /] for odd A/,. Such approx- 
imants are valid at any energy and are able to fit experimental 
phase shifts with high precision, even for a limited order. 

Expansion (|9j has the additional advantage that the corre- 
sponding single-channel potential is known analytically l6i]. 



Vi,MAr) = ^^-2—\r.W[u,, 



■,UMi\ , 



(11) 



where is a Wronskian determinant of functions Ujn{r)- 
These functions are solutions of the free radial Schrodin- 
ger equation, —u'^ + l{l + l)r~'^Um = EmU,n, at energies 
E„i = —i^lm. given by the iS-matrix poles. 

Combining this single-channel inversion technique with the 
EPP transformations provides a complete coupled-channel in- 
version scheme. First, experimental eigenphase shifts are in- 
verted separately for each channel, i.e., they are fitted with Eq. 
(|9|l and the corresponding potential is constructed with Eq. 
( fTTI ). This leads to diagonal potential and scattering matrices, 
which can then be used as a starting point for EPP transforma- 
tions. The factorization energies of these transformations are 
chosen by fitting the experimental mixing parameter with Eq. 
(O and the corresponding coupled-channel potential is finally 
constructed by iterated application of Eq. (O, followed by a 
compensation constant rotation. The final potential is thus di- 
rectly built from its scattering-matrix poles, which are either 
associated to the eigenchannels or to the coupling. 

Let us now consider one of the most important coupled- 
channel inverse problems, i.e., the inversion of the neutron- 
proton ^Si-^Di scattering matrix. To test our method, we 
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FIG. 1: (Color online) Poles of the neutron-proton '^S\-^D\ scatter- 
ing matrix close to the origin (note the different scale for the imagi- 
nary energy axis). 

use the scattering matrix of the Reid93 potential fl^ as in- 
put data. Using Eq. (|9]l with 5 terms for both the s and d 
waves {Ms — Md = 5), we find the following ten factor- 
ization constants s,„ = 0.23154, -0.45146, 0.43654, 1.6818, 
2.3106 fm-\ dm = -0.36719, -0.54420, 0.34828, 0.71766, 
3.3758 fm^^, which fit the neutron-proton ^Si and '^Di phase 
shifts, respectively. Note that si = 0.23154 fm~^ is fixed to 
reproduce the deuteron binding energy Bd = 2.2245 MeV. 
The S'-matrix poles — and — closest to the origin are 
represented in the complex energy plane in Fig. [T] 

In Fig. 121 we compare the eigenphase shifts of the Reid93 
potential and our fit. The coiTesponding potentials T4;5 and 
Vd;5 are determined by Eq. (fTTl i with the following sets 
of solutions: Us-m = exp(smr), Us-n = sinh(s„r), and 
Ud-m = [Sdmr cosh{dmr) - (3 + d^r^) sinh(dmr)] /r^, 
Ud;n = e-'*"'- [l + 3/(d„r) + 3/(d„r)2], m = 1,2, n = 
3, 4, 5. After extracting the centrifugal term from the second 
channel, we plot these potentials in Fig. [3] as well as their 
asymptotic behavior in logarithmic scale in Fig. |4] At large 
distances, the one-pion-exchange (OPE) short-range behavior 
of the Reid93 potential, with a pion mass of ^ 0.684 
fm~^, is clearly seen. Both inversion potentials are also short 
ranged. For the s-wave potential, this is surprising at first 
sight: in general, for a SUSYQM inversion, the leading con- 
tribution to the potential asymptotics comes from the pole 
iKm with Km > closcst to the origin, and behaves like 
exp(— 2Kmr). In the present case, this is the bound-state pole 
isi, which lies closer to the origin than the OPE cut, and one 
would expect the inversion potential to decrease slower than 
the OPE potential at large distances. We solve this problem 
here by exploiting the degree of freedom provided by the pres- 
ence of a bound state: in our inversion technique, the bound- 
state asymptotic normalization constant (ANC) can be cho- 
sen arbitrarily without affecting neither the phase shift nor the 
binding energy |11]. For a particular value of the ANC, re- 
lated to the residue of the S'-matrix bound-state pole lfl2ll . the 
exp(— 2sir) potential tail vanishes; for our single-channel S'- 
matrix parametrization, and thus for potential Vs-5, this ANC 



TABLE I: Real Erj and imaginary Eij parts of the P energies Ej 
entering the EPP transformations (j = 1, . . . , P). 

P,j 1,1 2,1 2,2 3,1 3,2 3,3 

-Erj (fm-2) 0.1463 0.1430 2.191 0.1590 4.004 
Eij (10-^ fm-2) 6.372 5.864 49.85 8.153 -3107 9443 



TABLE II: Deuteron properties (s-wave ANC, d-wave probability) 
for the Reid93 and the P-transformation inversion potentials. 

Property Reid93 P=l P = 2 P = 3 

As (fm-i/2) 0.8853 0.8851 0.8851 0.8851 
Pd (%) 5.699 5.27 5.39 5.73 



has the value As-5 = 0.8854 fm^^/^. 

Despite this satisfactory short-range behavior, potentials 
Vs-5 and Vd:5 differ from the Reid93 diagonal potentials: Figs. 
[3]and|4]show that Vs-5 is more attractive, whereas Vd-5 is more 
repulsive. This confirms that a coupling is necessary to re- 
cover the original potential. We thus transform the diagonal 
potential Vq = diag(ys;5, Vd-5) into a coupled one by using 
EPP transformations. The P complex energies Ej are cho- 
sen to reproduce physical properties of the Reid93 potential. 
For positive energies, we fit the mixing parameter, while for 
the negative bound-state energy, we fit the ratio of ANCs for 
the d- and s-wave components, r] = Ad /As, to its value for 
the Reid93 potential, 77 = 0.0251 |[T3l. These ANCs are re- 
lated to the residue of the two-channel scattering matrix at the 
bound-state pole [13] and read As ~ As-.5 cos e{—Bd), Ad — 
—As-.5sme{—Bd)- Hence, one has to satisfy tane(— i?^) = 
—0.0251. For the Reid93 potential, the mixing parameter is 
not an odd function of the energy, which confirms the need for 
the general EPP transformations introduced above. A single 
transformation (P = 1) provides two free parameters, Efn 
and En; we fix them to reproduce 77 and the slope of the mix- 
ing parameter at zero energy, e'(0) = 0.297 MeV^^. The val- 
ues of Ej^i, Eji compatible with these constraints are given in 
TableUfor P = 1. The mixing parameter £2 is shown in Fig.|2] 
by dots. It is seen that this model provides a good low-energy 
fit, in the spirit of the Newton-Fulton potential |2]; however, 
the present potential is much simpler and more general, as it 
allows a mixing parameter which is not an odd function of the 
energy. 

Fitting the mixing parameter at higher energies requires ad- 
ditional free parameters. Therefore, we iterate the EPP trans- 
formations. In Table H] we list the positions of the S-matrix 
poles for P = 1; 2, 3 EPP transformations. The correspond- 
ing mixing parameters are plotted in Fig. |2l which shows a 
good improvement with increasing P. The effective poten- 
tials are shown in Fig.[3](the P = 1 and P — 2 potentials are 
hardly distinguishable at the scale of the figure). They still dis- 
play some differences with the Reid93 potential, due to their 
different high-energy S-matrix, but Fig.|4]shows that, asymp- 
totically, our P = 3 potential is very close to the Reid93 po- 
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FIG. 2: Eigenphase shifts S3 g^, Ss^.^ and mixing parameter ei of the S'l- Di channels for neutron-proton scattering. 
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FIG. 3: Reid93 potential and its inversion potentials for the neutron-proton ''S'l-^Di channels. 
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FIG. 4: Asymptotic behavior of the Reid93 and inversion potentials. 



tential. In Fig.[T] we show that some of the S'-matrix poles are 
concentrated near E = — and are nearly independent 
of P. These poles can be associated with the OPE contribu- 
tion to the np interaction. Other deuteron properties are given 
in Table HIl for the inversion potentials; a quick convergence 
towards the Reid93 values is also seen. 

In conclusion, we have presented in this Letter an op- 
timal coupled-channel inversion algorithm, combining stan- 
dard single-channel SUSYQM transformations with coupled- 
channel EPP transformations that display the remarkable fea- 
ture of modifying the mixing parameter without affecting the 
eigenphase shifts. The method is developed for two chan- 
nels with equal thresholds and we plan to study its applica- 
bility to more general cases. In the neutron-proton ^Si-^Di 
case, the method leads to a simple potential, directly related to 



the position of its scattering-matrix poles. This very efficient 
parametrization could also be used in a phase-shift analysis of 
experimental data, directly providing the corresponding po- 
tential. 
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